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We consider the modification of the Cahn-Hilliard equation when a time delay process through 
a memory function is taken into account. We then study the process of spinodal decomposition in 
fast phase transitions associated with a conserved order parameter. The introduced memory effect 
plays an important role to obtain a finite group velocity. Then, we discuss the constraint for the 
parameters to satisfy causality. The memory effect is seen to affect the dynamics of phase transition 
at short times and has the effect of delaying, in a significant way, the process of rapid growth of the 
order parameter that follows a quench into the spinodal region. 
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I. BACKGROUND 

Diffusion is a typical relaxation process and appears in 
various fields of physics: thermal diffusion processes, spin 
diffusion processes, Brownian motions and so on. It is 
empirically known that the dynamics of these processes is 
approximately given by the diffusion equation. Although 
the diffusion equation has broad applicability, there exist 
the applicability limitations. First, the diffusion equation 
does not obey causality [l|. Let us consider the following 
telegraph equation, 



d 2 d 
T—m{x,t) + — n(x,t) 



dt 2 



dt 



DV 2 n(x,t) = 0. 



(1) 



Then, the propagation speed is given by v = ^/D/t. 
Thus, the propagation speed of the diffusion equation is 
infinite because the telegraph equation is reduced to the 
diffusion equation in the limit of vanishing r. Second, 
the diffusion equation does not satisfy exact relations, for 
instance, the Kramers-Kronig relation and the f-sum rule 
@, Q . By solving the Heisenberg equation of motion, the 
exact Laplace-Fourier transform of the time-evolution of 
a conserved number density is given by 

Sn LF (k, z) = -C(k) + l--( n (0)) eq + 0(1/ z 4 ), (2) 
z z a m 

where S means the fluctuations from the equilibrium 
value and C(k) represents the Fourier transform of the 
correlation function of the number density. Because of 
the Kramers-Kronig relation and the f-sum rule, the term 
proportional to 1/z 2 disappears and the coefficient of the 
second term is proportional to the equilibrium expecta- 
tion value of the total number (n(0)) eq . These are not 
satisfied if the coarse-grained dynamics of the number 
density is assumed to be given by the diffusion equation 



It is known that the problem of causality and the sum 
rules can be solved by using the telegraph equation in- 
stead of the diffusion equation. Interestingly enough, it 
has been shown recently that the coarse-grained equation 
derived by employing systematic coarse-grainings from 
the Heisenberg equation of motion is not the diffusion 
equation but the telegraph equation [H, 0] . 

The discussion so far is applicable to conserved quan- 
tities because the diffusion equation is a coarse-grained 
equation of conserved quantities. On the other hand, the 
corresponding equation for a non-conserved quantity is 
the time-dependent Ginzburg-Landau (TDGL) equation 
in the sense that the TDGL equation is a overdamping 
equation. The microscopic calculation [f| again shows 
that the relaxation phenomenon is accompanied by os- 
cillation and cannot be described by the TDGL equation 
as is the case with conserved quantities. As a matter 
of fact, the equation has a similar form to the telegraph 
equation. 

The effect discussed above is, in particular, important 
for fast processes where the time scale of memory is not 
short enough, because the telegraph equation can be ob- 
tained by introducing the memory effect to the diffusion 
equation. Such fast phase transitions are expected to 
have happened in early universe and most certainly also 
characterize the phase transitions expected to occur in 
the highly excited matter created in relativistic heavy- 
ion collisions. In the early universe, such situations may 
have happened when the typical microscopic time scales 
for relaxation, given by the inverse of the decay width as- 
sociated with particle dynamics, is larger than the Hub- 
ble time. This is a situation likely to be expected when 
describing GUT phase transitions or even the inflation- 
ary dynamics @. In relativistic heavy-ion collisions, one 
expected to learn about the QCD phase transition. For 
instance, from the hydrodynamic analysis of the freeze- 
out temperature in the most central Au-Au collisions at 
130 A GeV, the typical reaction time is given by around 
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10-20 fm Q. However, the characteristic time scale of 
the memory function in the Langevin equation which de- 
scribes the dynamics of the order parameter of the chiral 
phase transition is predicted to be about 1 fm, which is 
not enough short to be ignored [f| . 



II. MEMORY EFFECT IN SPINODAL 
DECOMPOSITION 

So far, we have discussed the difficulties with causality 
and memory effects, and sum rules associated with linear 
processes. Next, we consider such problems associated 
with a nonlinear process, known as spinodal decomposi- 
tion^) @. 

We consider the general Ginzburg-Landau (GL) Free 
energy, 



F(cj>) 



W) 2 



(3) 



where is a conserved order parameter. To describe 
the phase transition, the parameter b is proportional to 
T c — T with T c being the critical temperature of the as- 
sociated second order phase transition. The dynamics 
of conserved order parameters can be described by the 
Cahn-Hilliard (CH) equation. In the ordinary CH equa- 
tion, the irreversible current induced by the GL free en- 
ergy is given by 



3(x,t) 



-rv 



6F(<t>) 



(4) 



where T denotes a kind of Onsager coefficient. To take 
the memory effect into account, the current is generalized 
as follows, 



3(x,t) 



[ dsd 3 x'M(x - x', t - s)V x , 4irT 
ollowing 



s) 



(5) 



For simplicity, we assume the following memory function, 

r 



Mix, t) = —e 

7 



(6) 



Here, the typical memory time is characterized by 7. 
Substituting the current to the equation of continuity, 
we obtain the modified CH equation, 



' dt 2 ' 



7^r#M) + i-Mx,t) = rv 2 ^. (7) 



d_ 

dt 



At the late stage of the spinodal decomposition, the be- 
havior of the process is described by the linearized equa- 
tion around equilibrium order parameter <pQ = •Jb/a, 

IgpMk, t) + — c (fc, t) = -Tk 2 {ak 2 + 2b)Mk, i), (8) 

whose solution is tfi c (k,t) = Aj~ exp(\' + t) + Bfc exp(A'_i), 
where Ak and Bk are arbitrary constants and 



-1± >/l - 4 7 rfc 2 (afc 2 + 2b) 
27 ' 



(9) 



On the other hand, the solution for the correspond- 
ing ordinary CH equation (7 = 0) is 4> n c{k,t) = 
Ck exp[— Tfc 2 (26 + ak 2 )t\. One can easily see that there 
exists a critical momentum for the solution of the modi- 
fied CH equation, k' 2 c = [~b + ^Jb 2 + a/(4 7 r)]/a. Below 
the critical momentum, the solution shows overdamping 
behavior similar to the ordinary CH solution, while above 
the critical momentum, the damping is accompanied by 
the oscillatory fluctuation mode. Thus, we can consider 
that the critical momentum is a kind of a ultraviolet cut- 
off because the higher momentum modes have rapid os- 
cillations and cancel in computing averages. The prop- 
agation speed is characterized by the group velocity. In 
the modified CH equation, this is given by 



v{k) 



VtW+2) 



(10) 



where D = 2bT and £ = \fajb. At k = 0, the group 
velocity reads i>(0) = yU/7. The maximum group ve- 
locity is a monotonically increase function of k. Because 
of the momentum cutoff discussed above, the maximum 
group velocity is given by 



2D 1 + £ 2 /(2 7j D) 



7 i+^/i+e/WfD) 



(11) 



For causality, Eq. pip should be less than one. This 
leads to the constraint 7 > D. In other words, we can 
always find allowed values of parameters for which the 
spinodal decomposition is causal, under the constraint. 
In Fig. [TJ the group velocity is plotted as function of D 
for a fixed value £ = 1. The three lines corresponds to 
7 = D, D + 0.1 and D + 1, respectively. One can see that 
the group velocity is larger than one for any D at 7 = D. 
At 7 = D + 0.1, one can find D where the group velocity 
is less than one. At 7 = D + 1, the group velocity is less 
than one for any D. The modified CH equation satisfies 
causality in this sense. 

Let us see how this also applies for instance to the 
spinodal modes and then use these results e.g. for the 
problem of searching a signal of a phase transition in 
RHIC. By considering the fastest- growing mode of the 
SD, which is defined by 



d -1 + A 



dk 2 7 



= 0, 



we derive the time scale of the fastest mode to be 

20 



T c (k r ) 



(12) 



(13) 



When 7 is very small, the time scale is reduced to 
Tnc(k r ) = 8£ 2 /D, which agrees with the time scale of the 
SD without the effect of memory Q. Since the modes 
that give the result of Eq. (Ti~3")) are in fact the dominant 
spinodal modes and r„ c < r c , this reflects itself in an 
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FIG. 1: The group velocity as function of D for a fixed 
value £ = 1. The three lines are corresponding to 7 = D, 
D + 0.1 and D + 1, respectively. 
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FIG. 2: Volume average of <j>(i) as a function of dimensionless 
time i for the noncausal (solid) and causal (dashed) solutions. 
The two curves that exit the top of the plot are results from 
the linearized theory. 



overall delay of the time formation of domains, as de- 
scribed by the causal CH equation compared to the non- 
causal one, as the phase transition proceeds. This feature 
is confirmed by our simulations shown below. For a prob- 
lem like RHIC and a possible signature of a phase transi- 
tion coming from it, this difference in time scales can be 
very pronounced and lead to a striking effect that a sig- 
nal, like charge fluctuations and domain formation, can 
be so much delayed that possibly could not be observed 
in the current experiments. For instance, in RHIC, the 
correlation length is typically £ ~ 1 fm. Using also the 
relation between the parameters D and 7 obtained for a 
quark plasma, 7 = 31?, which is consistent with our con- 
straint condition obtained from Eq. (jlip , and considering 
D ~ 3.7 fm [l,[l0j], we obtain from Eq. flT3J) that r c ~ 6.1 
fm, which is to be compared with the result r nc ~ 2.2 
fm. This represents almost a 200% difference for the 
time scales for the starting of the growth of fluctuations 
in Eq. ([7]) as compared to the ordinary CH equation (for 
7 - 0). 

We have solved Eq. ([7]) numerically on a discrete 
spatial square lattice using a semi-implicit scheme in 
time, with a fast Fourier transform in the spatial coor- 
dinates @|. We have checked the stability of the results 
by changing lattice spacings and time steps. In addition, 
for 7 > 1 we have also used a leap-frog algorithm and 
the results obtained with both methods agreed very well. 
For the noncausal 7 = equation, we used as initial con- 
dition </>(x, t = 0) a random distribution in space with 
zero average and amplitude 10~ 3 . For 7 7^ 0, we used 
in addition the condition that at t = the first-order 
derivative of <f> is zero. For the numerical work, Eq. ([7]) is 
re-parameterized to dimensionless variables, conveniently 
defined by time t = (8/r„ c )t = (2b 2 T/a)t, space coordi- 
nates Xi — Xi/£, field (f) — \J ~cjb(j> and 7 = (8/r nc )7. 
In terms of these variables Eq. (O becomes function of 
only one parameter, 7. Eq. ((?]) was next solved for sev- 
eral values of 7. Two representatives results, for 7 = 



and 7 = 40.5 (for the example analyzed in the previous 
paragraph) , are shown in Fig. [2] 

In Fig. [2] we present the time evolution of the vol- 
ume average of the order parameter cf>, defined as <f>{t) — 
~W Sx(^( x ! *)+)> WIiere N 3 is the total number of lat- 
tice points and the average is taken over different initial 
random configurations. </>(x, i) + indicates that only the 
positive values of the field are considered (i.e., a specific 
direction for the field has been selected). The rapid in- 
crease of (f> reflects the phenomenon of SD. The figure also 
shows the results by solving the linear equation. It shows 
that it performs extremely well up to and right after the 
spinodal growth of the order parameter, then justifying 
our previous analytical results based on the solution of 
the linear equation. The effect of 7 is seen to be more 
important at earlier times, consistent with the memory 
function used and becomes less important after the rapid 
growth of the order parameter (the spinodal explosion). 
It also shows the effect of a finite 7, increasing dramat- 
ically the delay of the spinodal explosion, as predicted 
by our previous analytical results. The time for reaching 
equilibrium is seen to be very long, as is common with the 
traditional noncausal CH equation. Also apparent from 
Fig. [2] are the oscillations in the order parameter for a 
finite 7, also predicted by our previous analysis. This 
is due to the increasing importance of the second-order 
time derivative as compared to the first-order one as 7 in- 
creases, i.e. as 7 increases the dissipation term becomes 
less important and the equation becomes more and more 
a wave-like equation. The estimated delay for the ther- 
malization is even larger than the recent estimation [TlT | 
for the time delay of the relaxation of a nonconserved 
order parameter. 

We have also investigated the effect of memory for the 
structure function, S(k,t) = \4>(k,t)\ 2 . This quantity is 
important because it provides information on the space- 
time coarsening of the domains of the different phases. In 
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t 

FIG. 3: Spherically averaged structure function as a function 
of dimensionless time i for a k — 0.7k c for the noncausal 
(solid) and causal (dashed) solutions. 



be enough time for the system to thermalize before the 
breakdown of the system due to expansion. 

In this paper, we simply discussed the case of mem- 
ory function which has the exponential form. However, 
we can generalize the present work by considering other 
types of memory functions. In particular, microscopic 
calculations from nonequilibrium quantum field theory 
Q show that memory functions tend to exhibit relax- 
ation as well oscillations as time goes on. Then, we can 
consider the following physically motivated from those 
field theory calculations, 

M(x,t) = -e-* /7 cos(m),5 (3) (x), (14) 

7 

where fi is some characteristic frequency of oscillation. 
Substituting Q14p into the equation of continuity, we ob- 
tain 



Fig. [3]we present the results for the spherically averaged 
value of S for k — 0.7 k c (to emphasize the fast growth 
of the long wavelength fluctuations with k < k c ). The 
spherical averaging for a given k was done over momenta 

kr — 



k' 2 



U2 

Ky 



k 2 such that k - 0.1A < k r < k 



0.1 A, with A = 2n/L, where L is the size of the lattice. 
Consistently with Fig. [2j this figure shows the dramatic 
delay for the spinodal growth. 



III. CONCLUSION 

We have introduced memory effects into the CH equa- 
tion. By introducing the memory effects, we could de- 
fine the group velocity and, to satisfy causality, we were 
able derive the constraint for the parameters of the CH 
equation. For a physical situation of interest for the 
phenomenology of RHIC, we found that the inclusion 
of memory effects can delay substantially the phase- 
separation process and consequently, there might not 



(9 3 d 2 d 

7 2 ^(x, t) + 2 7 ^0(x, *) + (! + l 2 n 2 )-^ t) 



rv 



s<t> 



d_ 

at 



(15) 



This result and the example worked out explicitly in this 
work shows that the macroscopic dynamics of conserved 
quantities strongly depends on the choice of the memory 
function. In these cases, the use of microscopic motivated 
memory functions based on the particular model under 
study is fundamental. An analysis for (JTSJ), analogous to 
the one performed in this work for the simplest exponen- 
tial decay memory function, will be presented elsewhere. 



Acknowledgments 

The authors would like to thank CNPq, FAPERJ and 
FAPESP for the financial support. 



[1] D. Jou, J. Casas- Vazquez and G. Leben, Rep. Prog. 

Phys. 51, 1105 (1988); ibid. 62, 1035 (1999). 
[2] L. P. Kadanoff and P. C. Martin, Ann. Phys. (N.Y.) 24, 

419 (1963). 

[3] T. Koide, Phys. Rev. E72, 026135 (2005). 

[4] This microscopic calculation further makes it appear that 
the relaxation time r contained in the telegraph equation 
is not constant and has momentum dependence. 

[5] T. Koide and M. Maruyama, Nucl. Phys. A742, 95 
(2004). 

[6] M. Gleiser and R. O. Ramos, Phys. Rev. D50, 2441 



(1994); A. Berera and R. O. Ramos, Phys. Rev. D63, 

103509 (2001); Phys. Lett. B607, 1 (2005); Phys. Rev. 

D71, 023513 (2005). 
[7] Y. Hama, T. Kodama and O. Socolowski Jr., Braz. J. 

Phys. 35, 24 (2005). 
[8] T. Koide, G. Krein and R. O. Ramos, Phys. Lett. B 636, 

96 (2006). 

[9] D. Bower and S. Gavin, Phys. Rev. C 64, 051902 (2001) 
[10] M. A. Aziz and S. Gavin, Phys. Rev. C 70, 034905 (2004). 
[11] E. S. Fraga and G. Krein, Phys. Lett. B614, 181 (2005). 



